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Abstract 

Anomalous events with a large transverse momenta Pt were found 
for all secondary particles in an analysis of experimental data on cos- 
mic ray interactions at E'o > 5 • IQ^^eV . In order to clarify the under- 
lying mechanism of these events, the intermittency analysis of interac- 
tions with both large and usual transverse momentum was performed. 
It is shown that behavior of intermittency indices '^{q) as a function 
of rank factorial momentum q for such interactions is very peculiar. 
The indices show slight growth for the events with large transverse 
momenta for all secondary particles in contrast with the fast growth 
for events with normal average transverse momenta for final particles. 
The significant difference in the behavior of intermittency indices gives 
strong evidence on the existence of the new class of hadron interac- 
tions with large Pt for all secondary particles in the event. A possible 
interpretation of these anomalous events as a manifistation of Color 
Glass Condensate in high energy hadron-hadron interactions is given. 
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1 Introduction 



Investigation of the properties of quark-gluon interactions at TEV region of 
energy is very important task, especially in the light of future experiments 
at LHC One of the important goals of these experiments is clarification of 
the structure of high density/temperature quark-gluon matter which can be 
produced in heavy ion collisions. At present, cosmic ray experiments give 
the unique information on possible interesting effects in strong interaction of 
hadrons, which could be carefully investigated at LHC. 

A possible manifistation of Color Glass Condensate (CGC) in high energy 
collisions is now widely under discussion (see recent reviews P]). Phenomena 
of CGC arise in the high parton density regime of QCD where the typical 
virtualities of quarks and gluons are not very large and thus their interaction 
is not very strong. One of the most promising microscopical candidates 
for CGC is the sphaleron state ^2] which is closely related to the instanton 
model of QCD vacuum ( see reviews [H]) The feature of the sphaleron driven 
multiparticle production is rather high multiplicity in the events and large 
transverse momentum for all secondary particle. 

In this Letter we will discuss a possible evidence for the formation of 
CGC in high energy {Eq > 5 ■ cosmic ray interactions with nuclei 

in the stratosphere. An analysis of behavior of intermettency indices for un- 
usual type of the events, where all the secondary particles were produced 
with a large transverse momenta Pt ~ IGeV/c 4j is performed. An in- 
teresting feature of these events is a large difference in Pt distribution of 
the secondary particles compared to their distribution in the "usual" events, 
which have Pt ~ 0.2GeV/c. Similar events were observed by JACCE 
and Concorde [HI Collaborations in emulsion chambers. In spite of the fact 
that the first data on existence of "anomalous" events were obtained about 
twenty years ago, the underlying mechanism of hadron-hadron interaction, 
which is responsible for such events, is not understood so far. A promising 
way in investigation of dynamics of multiparticle production is the study of 
structure of the density fluctuations in the events through calculation of so- 
called intermettency indices (7|. The values of the indices are very sensitive 
to the microscopical mechanism of the hadron final state formation and, as 
we will show below, are very different for "anomalous" and "usual" cosmic 
rays events. 
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2 Experimental data 



In this Section we present a new analysis of the experimental data, obtained 
in cosmic ray interactions at Eq > 5 ■ which were observed in X-ray 

film - emulsion chambers exposed in the stratosphere at elevations > 25km. 
Each chamber contained target and gamma blocks. The target blocks were 
gathered of many layers of heavy or light material interplead with thin nuclear 
emulsions; gamma blocks contained 4-10 lead plates, 5mm thickness each 
with nuclear emulsions, and X-ray films between them. The details of the 
construction of the experimental chambers are described in 

For each event in the target block we determined the number of charged 
secondary particles (rich), the number of secondary 7 - quanta (n^), the angles 
of emission of each secondary particle in the laboratory frame, Qch, and, B-^ 
, an energy and a transverse momentum of each over threshold 7 - quantum 

and Pt-j (the threshold for the electromagnetic cascade in X-ray films was 
2TeV and in nuclear emulsion - O.OSTeV^). We estimate the deviation of the 
accuracy for P-r-y not greater than 30%. The number of secondary charged 
particles and their angular distribution provide an estimation of the energy 
of the primary particles which produce interactions in the target block of 
the chamber. The analysis of measured values of transverse momentum of 
each photon in events with > lOTeV showed that 7 of them (index 

" Pt = 1") have rather different properties in compare with an others. Thus, 
the transverse momenta of most photons in these 7 events were several times 
larger than the average transverse momentum of secondary photons in the 
"usual" type events with < Ptj >~ 0.2GeV/c (index "Pt = 0"). 

The integral distribution of transverse momentum of all the secondary 
photons is given in Fig.l and can be described by the superposition of two 
exponents: 

X^(> pt^) = Aiexp{-pt-y/poi) + A2exp{-pt^/po2), (1) 

where poi ~ 0.2GeV/c for the "usual" events and po2 ~ 0.8GeV/c for the 
"anomalous" ones. The main feature of these 7 "anomalous" events is that 
most of the secondary photons have transverse momentum pt^ > 0.5GeV/c. 
The difference in characteristics between the two classes of events in Fig.l is 
large large compared to the errors in the estimation of E^ and can not be 
explained by loss of the underthreshold photons and statistical fluctuations. 
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3 Analysis of the pseudo-rapidity distribu- 
tions of the secondary charged particles 



At first glance, both types of events have a similar pattern. Thus, a com- 
parison of x^/NDF fit with Gaussian distribution for events with higher 
mutiplicity is given in table 1. The histograms of the pseudo-rapidity distri- 
butions of the secondary charged particles for these events are presented on 
Fig. 2. All events show approximately similar behavior. 



Table 1. Multiplicities and x^/^DF of pseudo-rapidity 
distribution for large and small Pt events 





small Pt events 


largePr events 


N 


NCH 


x' 


NCH 


x' 


1 


290 


1.99 


675 


3.09 


2 


239 


2.19 


432 


1.32 


3 


175 


1.24 


235 


2.42 


4 


174 


3.16 


182 


2.25 



4 Analysis of the structure of correlations 

The analysis of the scaled factorial moments of the multiplicity distribution to 
isolate the dynamical density fiuctuations from the background of statistical 
fluctuations was proposed in Ref. [Hj: 

< nt^l > 

<F,{5)> = (2) 

where < n^'^^ >=< nm{nm — 1) . . . {rim — Q + ^) >■ 

In the presence of self-similar fluctuations of different sizes, the depen- 
dence of the moment < > on the size of the phase-space bin follows the 
power law: 

> =(l/5)^« (3) 

for 6^0. The positive constant ipq is the so-called intermittency exponent. 

The distribution (2) is not a smooth function, i.e. it contains sharp 
spikes and holes between particles in the phase-space. The observation of 
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the power law in a sufficiently large range of scales 5 indicates a self-similar 
fractal structure of the short range particle density fluctuations. 

In the opposite case of a smooth distribution (probability density is con- 
tinuous), the factorial moments are 

< Fq> ~ Const. (4) 

It should be mentioned that the high order moments resolve the tail of the 
multiplicity distribution. Thus they are very sensitive to the density fluctua- 
tions at various scales of 5 used in the analysis. Therefore, the study of fluc- 
tuations at various values of 5, especially at small 5 and large q can improve 
our phenomenological understanding of multiparticle production processes. 

The particle multiplicity distribution is usually studied for a sequence of 
phase space domains 5 by consecutive subdivision of an initial region A in 
M equal sub domains: 

5 = A/M 

In every subdomain m (m = 1, M) Um is the multiplicity of that bin. 

In order to improve the statistical accuracy in the experimental estimation 
of factorial moments, F^'s of individual cells (1) are averaged over events and 
over M cells ("vertical analysis"). The moments averaged vertically {i.e. over 
events) can be defined double average 

where < >= tt— V] rim is the multiplicity in bin m. 

evt 

We use a modified method of vertically averaging jT2] in which moments 
are averaged at the starting point of the original region A location 

ivi ^ step step^^evt evt ^ ""m ^ 

where N^t^p is the number of small {step/ A <C 1) steps of the starting position 
of the original region A in the area of pionization. The pseudo-rapidity 
7] = —In tgQ/2 of the secondary particles charge is used as the basic variable. 
The initial value of A is 4.0 and M = 40. More detailed analysis has been 
performed with A = 3.5, 3.8, and M = 35, 38, respectively. 
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5 The results 



The results are presented in Table 2 and in Fig. 3-5. The values of In F^'s as 
a function of —In/ 61] for order g = 2 to g = 8, are plotted for small pt events 
in Fig. 3 and for events with large Pt in Fig.4. 

For small pt events in Fig. 3 the rise of In Fq when increasing —In 5r) 
takes place for all q. For large Pt events (Fig.4) at all q any significant 6r) 
dependence is absent. 

The slopes (pq in the region 0.1 < In 6ri < 1.0 indicated by the solid lines 
in the figures are given in Table 2. Values of are essentially larger for the 
events with small Pt compared to the large Pt events. 

The comparison of the slopes cpq for both types of events is shown in Fig.5. 
For small Pt events the data are consistent with the intermittent behaviour 
(2), while for "anomalous" events the values of the slopes and their changing 
as function of q are essentially smaller. 

Table 2. Slopes (pq fitted over 0.1 < 5rj < 1.0 for large and small Pt 
events 

small Pt events large Pt events 

q 



'fq fq 



2 


0.100 ± 0.004 


0.068 ± 0.005 


3 


0.260 ± 0.014 


0.095 ± 0.010 


4 


0.310 ± 0.027 


0.094 ± 0.016 


5 


0.51 ± 0.05 


0.08 ± 0.02 


6 


0.66 ± 0.06 


0.10 ± 0.03 


7 


0.77 ± 0.09 


0.11 ± 0.04 


8 


1.29 ± 0.11 


0.13 ± 0.06 



6 Possible explanation of anomalous events 

Self-similar particle-density fluctuations have been observed in various reac- 
tions [Zj. Thus, the approximate scaling in respect to decreasing of domain 
size was found. In spite of the strong experimental evidence, the effect is still 
rather far from understanding based on fundamental QCD theory. One of 
the possible contributions, which might give strong influence on behavior of 
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scaling moments, is coming from Bose-Einstein correlations. However, these 
correlation should not be related to some conventional static source 0|. 

There is a clear indication on Pt - dependence of intermittency in different 
processes ^01; QI]- An analysis for particles with Pt smaller or larger than 
0.15 4- Q.?>GeV/c in the same events reveals a strong sensitivity to the value 
of transverse momentum. The results show that the slopes (pq increase by a 
factor 2-^4 for value of cut Pt = O.lSGe^. A similar but not that strong 
effect is observed if the Pt cut is moved to O.SOGeV^/c. 

Our results for small Pt events confirm the power-law behavior given by 
Eq.2. On a contrary, for the large Pt events, expression (3) looks as a very 
promising candidate for the behavior of the intermittency exponents. 

The results on Pt distribution, presented in Fig. 1, and the results re- 
lated to the intermittency ( Fig. 3-5 and Tables 1, 2), strongly indicate 
the existence of a new class of events with large Pt for all secondary par- 
ticles. We should mention that intermittency dependence for the "usual" 
events might be understood as a result of perturbative QCD evolution of 
partonic cascade at high energy. In this case the independent hadronization 
of the quarks and gluons in partonic chain should lead to a significant inter- 
mittency patterns [7j. For the "anomalous" large Pt events, we found that 
the behavior of the factorial moments is very different compared to that of 
the "usual" events. From our point of view this is very important result. 
It seems that in this case the perturbative QCD can not describe the data 
and a new mechanism for secondary particle production should be involved. 
One of the possible candidates for this mechanism is the contribution from 
nonperturbative partonic interactions at high energies, related to the Color 
Glass Condensate formation at large parton density P] . In this case the pro- 
cess of particle production in very high energy nucleus-nucleus interactions 
induced by cosmic rays looks like simultaneous burst of all interacted regions 
in the overlapped cosmic and target nuclei. If the CGC is identified with the 
QCD sphaleron state j2], which describes the particle production induced 
by strong non-perturbative fluctuations of vacuum gluon fields, instantons 
[Hj, than the transverse momentum of produced particles should be roughly 
PTilarge) ~ 1/pc ~ Q.QGeV , where pc ~ 0.3 fm is the average size of instan- 
tons in QCD vacuum. This value of transverse momentum is much larger 
than the value of transverse momentum PT{small) ~ l/Rconj ~ 0.2GeV 
{Rconf ~ Ifm is the confinement size) of secondary particles in "usual" 
events. 
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7 Conclusion 



In summary, a strong indication of the existence of a new class of hadron- 
hadron interactions at high energies is found. A possible explanation of these 
interactions is a manifistation of the Color Glass Condensate formation in 
the cosmic ray induced events. It would be interesting to search for events 
with similar properties at current (HERA, RHIC, TEVATRON) and future 
(LHC, TESLA) coUiders. 
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FIGURE CAPTIONS 



Fig. 1. Pt distribution of secondary particles in two types of events. 

Fig. 2. Pseudo-rapidity distribution of events with small and large Pt- 

Fig. 3. In < Fg > as a function of —In Sr) for the group of small Pt events. 

Fig. 4. In < Fq > as a function of —In 6 for the group of large Pt events. 

Fig. 5. The comparison of the slopes (pq as function of the moment order q 
for two types of events. 
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